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i Abstract. We carry out a detailed analysis of a time dispersive dissipative 

1-" ' (TDD) string, using our recently developed conservative and Hamiltonian ex- 

yJ ' tensions of TDD systems. This analysis of the TDD string includes, in par- 

MH , ticular: (i) an explicit construction of its conservative Hamiltonian extension, 

VQ ' consisting of the physical string coupled to "hidden strings;" (ii) explicit for- 
mulas for energy and momentum densities in the extended system, providing 

*vj ' a transparent physical picture accounting precisely for the dispersion and dis- 

►^ , sipation; (iii) the eigenmodes for the extended string system, which provide 

an eigenmode expansion for solutions to the TDD wave equation for the TDD 

^"^ ' string. In particular, we find that in an eigenmode for the extended system 

^~^ ' the displacement of the physical string does not satisfy the formal eigenvalue 

_^ , problem, but rather an equation with a source term resulting from the exci- 

/— ^ 1 tation of the hidden strings. The obtained results provide a solid basis for 

v^ ' the rigorous treatment of the long standing problem of scattering by a TDD 

^»^ ' scatterer, illustrated here by the computation of scattering states for a string 

■"".^ , with dissipation restricted to a half line. 



s 



1. Introduction 



Q , The need for a Hamiltonian description of dissipative systems has long been 

known, having been emphasized by Morse and Feshbach [5, Ch 3.2] forty years ago. 
Recently we have introduced conservative [Ij and then Hamiltonian [51 13] theories 
^^ ■ of time dispersive and dissipative (TDD) systems addressing that need. The ex- 

tended Hamiltonian is constructed by coupling the given TDD system to a system 
of "hidden strings" in a canonical way so that it has a transparent interpretation 
as the system energy. It turns out that such a system of hidden strings is, in fact, 
a canonical heat bath as described in 4, Section 2] , [5] Section 2] . 

This work is intended as a self contained companion to the papers [U [2 ISIj 
illustrating the phenomena described there through explicit calculations for the 
typical example of a TDD string. 

To keep everything elementary, let us consider a scalar wave equation in one 
spatial dimension, such as might be used to describe wave propagation along a 
homogeneous string. In the absence of dispersion, the displacement (f){x,t) of the 
string at position a; G M and time t g M evolves by the ID wave equation 

(1.1) dU{x,t) - idl^{x,t) = /(x,t), x,ieR, 

where 

• 7 is the tension of the string. 
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• we have taken units in which the mass per unit length is 1. 

• f{x, t) is an external driving force per unit length, usually supposed bounded 
and compactly supported. 

For the most part, we consider (|l.ip with the string at rest at t — — cx), 

(1.2) lim 4>{x,t) = lim 9(0(2;, t) = 0. 

t — ^ — oo t — *oo 

Thus the solution (f)(x,t) is a function of the driving force f{x,t), indeed 

(1.3) (/.(x,t) = - / / f{y,t')dyAt' 

^"^ J ~oc J x-v(t-t') 

with V — 1/a/7 the speed of propagation on the string. 

In this paper we consider a modification of p.ip incorporating friction in the form 
of a dissipative term with time dispersion. Specifically, we consider the following 
equation 

(1.4) dtUHx,t) + l x{^,t-t')dt'Hx,t')dt'\-^dlcj,{x,t) = f{x,t), 

with xix, t) a given function, called the susceptibility, satisfying a power dissipation 
condition — (|1.23p below. 

The physical idea behind (|1.4p is as follows. The wave equation (jl.ip can be 
expressed in terms of the string momentum tt{x, t) = dt(t>{x, t) as 

(1.5) dtn{x,t) = sum of all forces at a: = f{x,t) + ^d'^4>{x,t) 

driving force string tensile force 

In the damped string p.4p . the basic relation dtir = "sum of all forces" still holds, 
but the simple relationship n — dt4> between the string momentum and velocity is 
replaced by the material relation 

(1.6) 7r(x,i) ^ dt(t){x,t)+ I x{x,t~t')dt'<t){x,t')<lt'. 



Eq. (|1.6p is supposed to be a phenomenological relation describing the interaction 
of the string with a surrounding medium, expressing the fact that some of the 
string momentum is absorbed by the medium and then partially retransmitted to 
the string with delay. 

The wave equation (jl.ip may be expressed as a Hamiltonian system. As de- 
scribed in [51 12| , the procedure of going from p.ip to (|1.4p is naturally understood 
in this context and indeed (|1.4p can he derived from a larger Hamiltonian system, 
with additional variables . The additional variables of this extension may be inter- 
preted as describing a "hidden string," with internal coordinate s, attached to each 
point X of the physical string (see figure [J) . These "hidden strings" are coupled to 
the physical string via a coupling function q{x, s), but do not interact directly with 
one anotherl^ Upon "integrating out" the hidden strings, the result is an effective 



The lack of interaction between hidden strings is not essential. Indeed the methods of [T][2l|3] 
allow to consider systems with "spatial dispersion" in the material relation 

7r(x,t) = dt4>{x,t)+ / x{x' ,T)dt4>{x-x' ,t-T)dTdx' . 
Jr Jo 

The extensions for such systems involve interaction between the hidden strings, but for simplicity 

we do not pursue this here. 
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Figure 1. The extended system consists of the physical string 
(thick/red) coupled to an independent "hidden string" (thin/blue) 
at every point of dissipation. This Hamiltonian system models 
exactly the dissipation and dispersion in the physical string via 
the exchange of energy with the hidden strings. 



dynamics of the form p.4p for the physical string, with susceptibility 

1 
2 



(1.7) 



Xix.r) 



i,{s,x) 



S + T 



(,(x, r) dr ds. 



Different choices of the coupling function c^{x, s) produce different susceptibilities 
x(a;,T), and furthermore [21 [^ any susceptibility satisfying the power dissipation 
condition — (|1.23p below — can be obtained in this way. 

Given the role of the TDD wave equation (|1.4p as a phenomenological description 
for the interaction of the string with its environment, it is not surprising that 
an extended system exists which produces the effective dynamics (|1.4I) . at least 
approximately. However, the main point of [51[3], which we explain below, is that 
it is possible to construct explicitly an extended system which exactly reproduces 
(|1.4[) . Further, using that system we can derive simple expressions for quantities 
like the energy and wave momentum densities for the damped string in terms of 
the susceptibility x- 

In addition, as illustrated below, the extended system clarifies the nature and 
role of eigenfunctions in a TDD system by providing a natural eigenmode expansion 
for solutions to (|1.4p . In this context, an eigenfunction exp(—iujt)(j)^{x) is a time 
periodic steady state solution to (|1.4p . For the simple non-dissipative string (|l.ip . 
an eigenmode solves 



(1.8) 



-idl(t)^{x) 



'Uj'^(t>^{x), 



which is to say it is a linear superposition of the plane waves exp(±ia;a;/^/7). For the 
TDD string ()1.4|) . the dissipation induced by the dispersion in ()1.6|) may preclude 
the existence of a steady state solution without a source term. However, there is no 
difficulty in constructing the eigenmodes for the extended system, as it is without 
dissipation. We find below that the resulting eigenmode equation for (f)^ 

(1.9) ^dl(j)^{x) = -uj"^ {l + Rex{x,uj))(j)^{x)~\ujHvax{x,uj)g^{x), 
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has a source term involving an arbitrary function g^^ , which depends on the config- 
uration of the hidden strings. Here x is the r- Fourier transform of %, 

/"OO 

(1.10) x(a^,c^) = / e'^^x{x,r)dT. 

Jo 

Let us now sketch the construction of the Hamiltonian extension. To start we 

write (|l.ip as a first order system, introducing the momentum density tt{x, t) — 

dt(j){x,t), so that 

dt(t>{x,t) = ■n{x,t) 

^^■^^^ 9t^(x,t) = ^dl,p{x,t) + f{x,t). 

This first order system is Hamiltonian, with symplectic form 



(1.12) ^((/'i,7ri;02,7r2) = j {cJ)i{x)tt2{x) - (J)2{x)tti{x)) dx 

JR 

and time dependent Hamilton function 

(1.13) Hfi<f>,Tr,t) ^ \ [ {^{xf+-f\d^c^{x)\^) dx~ f fix,t)(f>{x)dx. 

^ Jr Jr 

As the system is non- autonomous, i//(7r(-, i), (/)(•, i), i) is not conserved. Instead 
along a trajectory to (|l.lip . we have 

(1.14) ^ = dtHfW-,t),7r{;t),t) = - f dtfix,t)(f>{x,t)dx. 

We take the internal energy of the (undamped) string to be the Hamilton function 
Hq{(J),tt) of the system with / = 0. Thus, using (|1.14p we obtain 

(1.15) ^ - ^ + i /'/(a;,i)0(x,t)dx = f f{x,t)dtHx,t)dx. 

dt dt dt Jr Jr 

In other words, the rate of work done on the system is the integral of the external 
force times the velocity - "power — force x velocity." 
Observe that the Hamiltonian Hq can be written 

1 

2 
where 

(1.17) K^ = y/^d^, K^ = Identity. 

Based on this expression, we separate the equations of motion (jl.lip , into "material 
relations" 

(1.18) f^{x) = K^c^ix) = ^d^cj){x) 

(1.19) /^(x) == K^Ti{x) = 7r(a;), 

and dynamical equations, 

dtd^{x,t) = KlUix,t) 
1.20 , 

dMx,t) = ~KlU{x,t)+f{x,t). 

To obtain a first order TDD system leading to (|1.4p , we follow [3J [3] and replace 
the identity (|1.19p with the "material relation," 

/>oo 
(1.21) A(x,t)+ / x{x,T)U{x,t-T)dT = TT{x,t), 

Jo 



(1.16) ifo(</',^) = 77 ll^-'^lli^(M) + ll^^<^lli^( 
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but maintain ()1.18|) and the equations of motion ()1.20p . i.e., 

dtcj){x,t) = Uix,t) 
(1.22) 

dtTT{x,t) = ■yd^<j){x,t)+f{x,t). 



(One could also introduce dispersion in the relation between and f^, or even 
dispersion mixing /^ and /^r, but for simplicity we consider here only TDD in the 
momentum.) The key requirements on the susceptibility x are 

(1) Causality, manifest in (|1.2ip in that the integral on the l.h.s. depends only 
on the past. 

(2) Power dissipation, which is the requirement that for every x 



(1.23) / g{t) x{x,r)dtgit-T)dTdt > 

J-oo Jo 

for an arbitrary function g{t), say compactly supported. 

The significance of the power dissipation condition is as follows. We suppose 
that the internal energy of the TDD system at time t is given by 

(1.24) i/o = 1 /" {U{x,tf+jd^(t>ix,tf) dx, 

^ Jr 

with fr,{x,t) as in (fr2T|) . Then 



jdx(j){x,t)dtdx(j){x,t) + U{x,t)dtU{x,t)jdx 
(jda:4'{x,t)dxf7:{x,t) + f.^{x,t)dtn{x,t) 
(1.25) - U{x,t) x{x,T)dtU{x,t~T)dTJdx 

(jdx<j)ix,t)dxU{x,t) + jU{x,t)dl(j){x,t) 
+ U{x,t)f{x,t) ~ Uix,t) xix,T)dtUix,t~T)dTJ dx 



f{x,t)dt(t){x,t)dx 





2, 



dt<j}{x, t) dx, 



X{x, T)dt (j){x, t-T)dT 
'0 

where in going from the third to the final line we have used integration by parts and 
the assumption that 8x4' ^'^d /^ vanish at spatial infinity. The first term on the r.h.s. 
is the rate of work done on the system by the external force /. Similarly, we interpret 
the second term as the power loss of the dissipative force, — / x{T)df(j){x, t — T)dr. 
Thus, power dissipation amounts to the requirement that, for any trajectory, the 
total work done by the dissipative force is negative. 

Let us note that the right hand side of p.23p is equal to 



(1.26) x{x,0) I 9{tydt+ g{t) drX{x,T)g{t-T)dTdt 

J -oo Jo 

\x{x, 0)5{ti - <2) + ^[drX]{x, \ti - t2\)}9{ti)g{t2) dH, 
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which foUows from integration by parts after noting that dtg{t — t) = —drgit — r). 
Here S{t) is the Dirac deha function. Thus the power dissipation condition is 
equivalent to the statement that for every x the (generahzed) function 



(1.27) D^{x,t) = x{x,0)6{t) + -[drx]{x,\t\) 



^1 

is positive definite in the sense of the classical Bochner's Theorem, see J71 Theorem 
IX. 9]. Thus, by Bochner's Theorem, the time Fourier transform of D^. is a non- 
negative measure, 

(1.28) D^ix,u;) = xix,0) + - [drx]{x,\t\)e'^' dt 

^ J-oo 

/•OO 

= x(a;,0)+ / drx{x,T)cos{ujT)dT 



— uj x(x,t) sin(aJT) dr — ujIuix{x,uj) > 0. 
Jo 

Since x is real, D-^(x,ll!) is symmetric under uj <-^ —to, and indeed Bochner's The- 
orem shows that a symmetric measure of suitably bounded growth is non-negative 
if and only if it is the Fourier transform of a real positive definite distribution. 

The simplest physically relevant example of a susceptibility satisfying the power 
dissipation condition is obtained with x(x,t) = a > 0, a positive constant. Then 



(1.29) Uix,t)+a f^{x,t-T)dT = TT{x,t). 

Using the equation of motion f^ = dt4>, wc find that 

/•OO 

■K{x,t) = dt(j){x,t)+a j dt(j){x,t-T)dT 
(i-.oO) Jo 

— dt(t>{x,t) + a(j){x,t), 

where we have applied the boundary condition limt^_oo 0(a;, t) — 0. Combined 
with the equation of motion dtTT = jd'^cj) + f we obtain 

(1.31) df(j){x,t)+adt^{x,t)~-fdl(j){x,t) = f{x,t), 

which is the dynamical equation for a driven damped string, with damping force 
per unit length —adt(t>{x,t). Note that a is dimensionally an inverse time and 1/a 
is the characteristic time for the damping of oscillations. 

A more realistic model for damping is obtained by supposing xC^^, t) to be a 
non-trivial function of r as in (ll.2ip . To allow for damping restricted to only a part 
of the string, we suppose that x depends on x as well. For instance, we could take 
the Debye susceptibility 

(1.32) x{x,r) = a(a:)e-^(-)^ 

with a{x) > and v^x) > non-negative functions of x. This results in the 
following integral-differential equation for </> 

/>oo 

(1.33) df(f>{x,t)+a{x) e-'^''^''df(f>{x,t^T)dT~jd^(j>{x,t) = f{x,t), 

Jo 
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or after integration by parts 

(1.34) at^0(x, t) + a{x)dt(l){x, t) - a{x)v{x)(t){x, t) 

+ a{x)iy{xf / e-''(^)^0(x,i-r)dT-7(9^0(x,t) = f{x,t). 



Jo 
Then 

- 1"°° UJ^ 
(1.35) D^{x,uj) = a{x) / e'^'^^^^^'w sin(tJT) dr = a(x)—-—^ j > 0, 

so the Debye susceptibility satisfies the power dissipation condition. 

The results of [U [3] show that under the above conditions it is possible to find 
a coupling function c;^(x, s), s e M, such that solutions to the TDD equations are 
generated by solutions to the following extended system 

dt(t>{x,t) = frr{x,t) 

dtTT{x,t) = jdl(t>{x,t)+f{x,t) 

dfipix, s,t) = 6{x,s,t) 
dt9{x,s,t) = 9fV(a;, s,t) +<^(x,s)/^(a;,t), 
with 

/oo 
<^(x, s)tp{x^ s,t) ds. 
-oo 

That is, given a solution (</>, tt, tp, 9) to (ll.36|) . at rest at i = — oo and with /tt given 
by p.37p , the first two coordinates (0, tt) obey ()1.22p with /^ given by ()1.21|) , and 
conversely any solution to (|1.2H ll.22|) with the string at rest at i = — oo arises in 
this way. Note that for each fixed x the additional variables 0, ip may be interpreted 
as describing the oscillations of a "hidden string" with coordinate s, displacement 
'ip{x,s,t), momentum 9{x,s,t) and driven by external force q{x,s)f.,r{x,t). Thus 
(|1.36p is precisely the system of coupled strings described above and illustrated in 
Figure [T] 

The extended system, consisting of the physical and hidden strings is Hamilton- 
ian with symplcctic form 



(1.38) >7(0i,7ri,'0i,6'i;(/)2,7r2,'02,6'2) = / {(/)i(a;)7r2(a:) - (/)2(a;)7ri(a;)}da; 

JR 

+ / {il)i{x,s)92{x,s) — ilj2{x,s)9i{x,s)}dsdx 

Jm? 

and Hamilton function 

(1.39) Hf{<j>,7r,^,,9,t) = i f {U{xf+-f{d^(t){x)f}dx^ f f{x,t)<i>{x)Ax 

^ JR Jm 

+ 1 f {9{x,s)^ + (dsi'{x,s)f}dsdx, 

where /^ is given by (|1.37p . that is f-^ix) — 7r(x) — / <;(a;, s)ip{x, s) ds. 

The dynamical equation for an excitation of the hidden string at x is a driven 
wave equation 

(1.40) d^'tPix,s,t) = d^,i;{x,s,t)+^{x,s)U{x,t). 
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The solution to (|1.40|) with the hidden string at rest at i = — oo is easily written 
down, see ()1.3p : 

1 pCC PS + T 

(1.41) ij{x,s,t) = - / <;{x,r)drU{x,t-T)dT. 

^ Jo J S-T 

Thus 

(1.42) / ^{x,s)ii){x,s,t)<ls 



1 f r+^ 

- / <;(x, s) / <;(a;,r)drds 

2 JR Js-T 



/^(x,i-T)dT. 



'0 

Comparing (|1.21|) and (|1.37p we see that for the extension to reproduce the TDD 
system upon elimination of ■0 and it is necessary and sufficient that the coupling 
<;(a;, s) satisfy (|1.7p . that is 



1 r /•^+''' 



(1.43) x{x,t) = - / <;{s,x) I (,{x,r)drds. 

The existence of such a function, which is unique under a natural symmetry as- 
sumption, is guaranteed by the power dissipation condition [21 [3] . 
Indeed, if we let ?(x, a) denote the s-Fourier transform of <j, 



(1.44) ?(x,f7) = / <?(x, s)e'"^ ds, 

then (ll.43|) is equivalent to 

(1.45) ^xi^^^) = ?(a;,^)?(2;,-^), 
so it suffices to take [2] 



(1.46) ^(x,g) = \J2D^{x,G). 

Furthermore, this choice of f(a;, a) is unique if we ask further that f(a;, cr) > and 
that a i—j ?(a;, cr) be symmetric. Then, 

(1.47) <;(x,s) = -^ f J2D^ix,a)e-'^'da = -^ f J2D^{x,a)cos{as)da, 

is real, symmetric, and positive definite. 

For example, in view of (|1.28p . (|1.35p . and (|1.47p . the following coupling produces 
the Debye susceptibility, 



(1.48) <;{x,s) = ^2a{x)ds'^{i^{x)s). 

where 



(1.49) *(s) = — / -\ -^e'^'dcu 



s I f^ 1 .,,,, . s 1 



ze 



"1*1 du = — - / e-l''l"'"("^) 



\s\TrJo sj\-v? |s|7r 

That is, 

(1.50) ^(x,s) = y2^(^(s)-^ / %in(0)e-''(")l^l™'('^) 
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2. Local energy and momentum conservation in the extended system 

We interpret the Hamiltonian Tin with / = as the internal energy of the 
damped string system consisting of the coupled physical and hidden strings. We 
have conservation of energy in the extended system, in the form 

(2.1) %Ho = f .fix,t)dtH^,t)dx, 

i.e., the rate of change of Ho is the rate of work done on the system by the external 
force. 

A significant advantage of working with the extended system is a transparent 
interpretation of the energy of the dissipative string as a sum of contributions from 
the physical and hidden strings. That is, it is natural to break the Tio into two 
pieces 

(2.2) noi(t^,7r,^,9,t) = Hoi^,7r,^)+Hhs{^,0), 
the energy of the physical string and hidden strings respectively, 

(2.3) i^o(0,^,V') = I f {U{x)^+j{d,(t){x)f}dx 

(2.4) Hhs{ij,e) ^\l {e{x,sf + {ds^p{x,s)f]dsdx, 

with /tt given by (|1.37p . 

The internal energy Tip can be written as the integral of a local energy density 

(2.5) £(x,i) = Eo{x,t)^Ehs{x,t) 
with 

Eo{x,t) = \{(dt<i>{x,t)f ^^{dA{xM) 

(2.6) { 

Ehsix,t) = - / {{dti^{x,s,t))^ + {ds^j{x,s,t))^}ds. 

^ JR 

The energy conservation law (|2.ip has the following local expression 

(2.7) dt£ix,t)+d,Jix,t) = f{x,t)dt4>{x,t) 
with the energy current 

(2.8) J{x,t) = --i,f^{x,t)d^<f){x,t) = --idt<f){x,t)d^<j){x,t). 

It is interesting to compute the time derivatives of Eo and Ehs alone: 

(2.9) dtEo{x,t)+d,J{x,t) = f{x,t)dt(t>{x,t)-dt/:^{x,t)dt(f){x,t) 

(2.10) dtEhs{x,t) = dt/^{x,t)dt(l){x,t) 
with 

(2.11) A(a;,t) = T:{x,t) ~ f^{x,t) 

/>oo 

<,{x,s)ip(x,s,t)ds = / x(2;,T)9t0(a;,i - r)dr, 

where we have used ()1.42p . The first of these equations p.9p is simply the local 
version of the energy law for the TDD string (|1.25p . From the second (|2.10p . we 
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see that the energy of the hidden strings, which is the energy lost to dissipation up 
to time t, is 

/t j*00 

dt'(j){x,t') / x{x,r)d^'<t>{x,t' ~T)dTdt'. 

= / / D-^{x,ti-t2)dt^(j){x,ti)dt2(j){x,t2)dtiAt2. 

J — oo J —oo 

If the susceptibihty x — ^^d hence the couphng <;^ — is independent of x, then 
the extended system is invariant under spatial translations. Associated to this 
symmetry is a local conservation law 

(2.13) dtv{x,t)^d^T{x,t) = -/(a;,i)9,0(x,t), 
for the wave momentum density 

p{x,t) = —7r{x, t)dx(t){x,t) — / 9{x,s,t)dxip{x,s,t)ds 

(2.14) -^^ 

= -{dt(j){x,t) + A{x,t)}dx(t){x,t)- / dtip{x,s,t)dxip{x,s,t)ds, 

JR 

with wave momentum flux, called stress, 

(2.15) T{x,t) ^ Ea{x,t) + A{x,t)dt(f>{x,t) 

+ j {{dt^{x,s,t)f^{ds4'{x,s,t)f]ds. 
When the driving force vanishes, / = 0, the total wave momentum 

(2.16) P = p{x,t)dx 



is a conserved quantity. 

The wave momentum density p ~ po + Phs is again a sum of contributions 

(2.17) pQ{x,t) = -{at0(x,t) + A(x,t)}a,0(x,t) 

(2.18) Phs{x,t) = - / dtil^{x,s,t)dxi^{x,s,t)ds 

from the physical and hidden strings. Likewise we separate the stress 

(2.19) T{x,t) = To{x,t)+ns{x,t) 
into two pieces, 

(2.20) To{x,t) ^ Eo{x,t) + A{x,t)dt(f>{x,t) 

(2.21) ns{x,t) - ^ f {idt^{x,s,t)f - idsyjix,s,t)f}ds. 

Then 

(2.22) dtPo{x,t)+dxTo{x,t) = -f{x,t)dx<p{x,t)+dxA{x,t)dt(f>{x,t) 

(2.23) dtPhs{x,t)+dxThs{x,t) = -dxA{x,t)dt(j){x,t). 

When we study eigcnfunctions below, it will be convenient to work with complex 
valued solutions. In the above expressions, terms which are quadratic in the field 
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variables should be modified in the complex case by the replacement ab -^ Re ab. 
That is 

1 



(2.24) 



etc. 



Eoix,t) = t{\d,^{x,t)\^+j\d,^ix,t)\^} 
Ehs{x,t) - ]- ( {\dti^{x,s,t)\'' + \dsi^{x,s,t)\'']ds 

^ .IB 



J{x,t) = —jRedt4>{x,t)dx4iix,t) 



Po{x,t) = -Redt(f>{x,t)da:4>ix,t) -ReA{x,t)dx:4>ix,t), 



(3.1) 



3. The eigenfunction equation 

It is useful and interesting to study steady state solutions to the extended sys- 
tem ()1.36p . for example solutions which are periodic in time e^"^*(0„(a;), 7r„(a:), 
il'uj{x,s),0uj{x,s)). We refer to the spatial component <I>i^(a;,s) = (0tj(a:), 7r^(a;), 
ipu{x, s), 6'(j(x, s)) of such a time periodic solution as an eigenfunction for the linear 
system (jl.36p with eigenvalue lo. Thus an eigenfunction satisfies 

-iuj4)^{x) = f^{x) 

-iuTT^ix) = ^dl(t>uj{x) 

-iujil;uj{x,s) = 9^(x,s) 

-iuj9uj{x,s) = d^ip^{x,s) +i;{x,s)U{x), 
with 

(3.2) 

We see that the displacement of the hidden string at position x satisfies 

(3.3) - d'^i>^{x,s) -uj^i>^{x,s) = ~iuj(,{x,s)(l)^{x), 
so 

(3.4) i;Ux,s) 

e-'"^^ ^?ix,a)da 

rr^ — - 1 .1^ 



fnix) + / <;{x,s)i)i^{x,s)ds - n^{x) = 0. 



a(x) cos[lljs) + b[x) sin(a;s) P.V. 



a{x) cos{lus) + b{x) sin(a;s) 



2tt 
T^(f>uj (x) 



sin(ci;|s' — s|)<;(x, s') ds', 



with a{x) and b{x) undetermined functions of x. Here, P.V. denotes the "principle 
value" integral, 



(3.5) P.V. 



1 



a^ 



LU^ 



y^{x, a)da 



= lim 

= lim — / e 

Sio 2 



^^-UJ^^\>S} 



j^{x, a) da 
1 



a^-{LU + iS)^ a-^-{u;~iSy 



^{x, a) da. 
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By (|3.4|) . we see that 
(3.6) / ^{x,s)ip^{x,s)ds = a{x)^x,uj 

JR 



27r 



/oo 
-^ T<r(a;, o')<;(x, -cr) da. 
-oo cr - ^ 



Recalling p.45p . that ^{x,a)'^{x,—a) — D^{x,a) = 2crlm x(a;, cr), we see that the 
final term on the r.h.s. can be expressed 



(3.7) -%(^P.V.r^^aImx(.,a)da 
_ i0(^(a;) 



P.V. / <^ + \alT^x{x.cT)da 

= -iujRex{x,uj)(f)^{x), 

where we have used the symmetry aim £(2;, cr) = — crlmx(a;, — cr) and the "Kramers- 
Kronig relation" 

1 /"^ 1 

(3.8) wRex(a;,w) = P.V.- / crlmx(a;, ct) dcr. 

I" J-00 cr-w 

Plugging this result into (|3.2p and combining the first two equations of (|3.ip we 
find that the string displacement (/) solves the following elliptic problem 

(3.9) jdlcjj^ix) = - u)'^ {1 + Rcxix, uj)) (f)^ (x) ~iLJa{x)tix,uj). 

In analyzing ()3.9p we should distinguish two cases: 

(1) f(a;, a;) = for all x. 

(2) f(a;,a;) non zero for x in an open set. 

In the first case, which is non-generic, the medium described by the hidden strings 
is not absorbing at the given frequency lu, that is x{xi'-^) is real for every x. Thus 
(13. 9p is a very strong restriction on 0^^ , namely that it should satisfy the Schrodinger 
equation 

(3.10) i-jdl + Vix))(bM) = A0„(x), 

with potential V{x) = — aj^x(a:,a;) and spectral parameter A = w^. In the second 
case, the physical string displacement may be decomposed as follows 

(3.11) M^) = (/.«(x) + (/.£2)(2;), 

where 0L {x) is any solution to the non-dissipative eigenfunction equation 

(3.12) ldl4,i'Hx) = ^u;\l + Rex{x,iu))^l'\x), 

and (f)uj is an arbitrary function with support in {x : 'q{x,uj) 7^ 0}. Indeed, given 
ipui and (pcj we need only choose a to be 

,,,., ,, . jd^Ju\x)+u;^{l+Rex{x,L^))cj,^J\x) 

(3.13) aix) = 1 ;;- 

UJ(;(X,UJ) 

to obtain a solution to (|3.9p . 

Thus, the physical string displacement for an eigenmode 4>ujix) may be chosen 
arbitrarily within that part of the string which is absorbing at the given frequency 
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Lo. We want to emphasize the significance of this fact, since the formal eigenvalue 
problem 

(3.14) ldlct>u.{x) - -uj\\+x{x,oj))ci>^{x) 

does not allow for such arbitrariness in the choice of (1)^1 {x) . The resolution to this 
apparent contradiction lies in recognizing that a TDD system is an open part of 
a larger conservative Hamiltonian system. It is only for the extended Hamiltonian 
system that the eigenmodes ^^ are unambiguously defined with e"^*$ij a stationary 
solution to the canonical Hamiltonian evolution equation. The "legitimate" eigen- 
modes for the original TDD string consist of projections (j)^ of the eigenmodes $1^ 
onto the subspace of the physical string. Thus a TDD string, being an open system 
has as many eigenmodes, i.e stationary solutions, as its the minimal conservative 
extension, introduced and described in [ll[2l[3]. In particular, a finite-dimensional 
TDD system typically has infinitely many stationary solutions and, hence, infinitely 
many eigenmodes. Another view on the construction of eigenmodes follows. 

The eigenfunctions written down above, involving as they do arbitrary excita- 
tions of the hidden strings, may not be relevant to the dynamics (ll.36|) with the 
external force acting on the physical string. Indeed, note that the effective equa- 
tion (|3.9p for (f)^ does not depend on the term b{x)sm{uos) appearing in (|3.4p . We 
shall see that eigenmodes with h{x) ^ are not needed in the expansion of a gen- 
eral solution to (|1.36p with a compactly supported external force. Essentially this 
is due to the fact that the configuration of the hidden strings remains symmetric 
throughout the evolution (|1.36[) . 

To proceed, it is convenient to introduce the Fourier-Laplace transform 

(3.15) h{C) = I e'^' hit) At, 

defined for 

(l)CGCif/i— >0 super exponentially fast as \t\ —> oo, for instance if h is 
compactly supported in time. 

(2) ImC, > ii h ^ super exponentially fast as t — ^ — cx). 

(3) ImC, < if /i ^ super exponentially fast as t ^ cx). 

Taking the Fourier-Laplace transform of (|1.36p results in the following equations 



(3.16) 



-iC(/>(x,C) = frrixX) 

-KnixX) = 7dl^ixX)+IixX) 

-iCV'(a;, s, C) = 0{x,s,O 

-iCe{x,s,C) = d^i;{x,s,C)+'i{x,s)U{x,C), 



with 

/oo 

<;(x, s)V'(x, s, C) ds. 
-oo 

Due to our consideration of solutions with the strings at rest at i = — oo, we 
expect the quantities in (j3.16|) to be well defined only for Im^ > 0. However, 
if the driving force is compactly supported in time then f{x, C) is defined for all 
C e C, and we may extend (f>{x, C) to C G C \ R by solving p.l6p . Following the 
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eigenfunction analysis, we obtain, with e = sign of IvciC,, 

27r lo, ^2 - C 



(3.18) V'(x,g,C) = - 'T / e-'-^:^-— ?(a:,a)da 



' — oo 



^(|C) /°°e'^';i^'-^l,(5',a;)d.'. 



(3.19) 75,'0(a;,C) = -C{^ + x{x,e(:))cj,{x,0 + .f{xX), 
where x(^7 C) i^ defined for Im^ > as 

/•oo 

(3.20) x(a^,C) = / e'^^x{x,T)dT. 

For a function h vanishing at +00 or — cx3 we have the Fourier inversion formula 

(3.21) h{t) ^ — / e-''^*=F'jt;^(w T i^) dw, f lim /i(f) = oV 

2tt J^ \t=±oo ) 

with ?7 > arbitrary. Inverting the solution to (|3.16p with ?7 > or < produces 
distinct solutions to (|1.36p : for ry > we obtain the desired causal solution with 
the strings at rest at t = —00, while for r; < we obtain the anti- causal solution 
with the strings at rest at i = +00. 

In a certain sense we are only interested in the causal solution to (|1.36p and 
thus to the solution to (|3.16p only for ( in the upper half plane. However, since 
p.l6p involves a source term, this solution, often called the resolvent, is not directly 
expressed as a superposition of eigenfunctions. However, there is a general proce- 
dure for decomposing the solution to (|3.16p into a superposition of eigenfunctions. 
Namely for w e M we define 

(3.22) ((/)^(x), 7r^(a;), V'u;(x, s), d^{x, s)) 

= lim ((j){x,uj + iS) — (t){x,uj — iS),n{x,uj + i5) — n{x,uj ~ iS), 

6 — *0 V 

ip{x, s,uj + iS) — ip{x, s,Lu — iS), 9{x, s,uj + iS) — cf){x, s,uj — i5)\ . 

Since f{x,() is continuous at each ^ = ^ G IR, it follows from (|3.16p that {(t)uj,Truj, 
V'uJi^uj) is an eigenfunction for each uj. To recover the resolvent for C in the up- 
per half plane from the eigenfunctions (|3.22p . suppose that the external force is 
supported in the set {i : t > to}- Then 

(3.23) |/(x,C)| < Ce-i-';*", (ImC>0). 

Thus exp(— iCio)/(a;, C) is bounded in the upper half plane, and by analyticity we 
have 13 

(3.24) (^(a:, C), ^(x, C), ^{x, s, C), 0{x, s, C)) 

1 f 0''^'^^'^''*° 

u{x),T:^{x),ij^{x,s),6^{x,s)) duj, {lm( > 0) 



27ri 7jj, uj-C 



For a general linear system the limit on the right hand side of 113.221 1 is a defined only as 
a vector valued measure and the integral on the right hand side of 1 13. 241 1 should be interpreted 
as the integral of l/{uj — Q against this measure. For the systems of coupled strings considered 
here, however, these boundary measures are always absolutely continuous, so 113.241 1 holds with 
{(/'u, TTij, V'lJi^ij) defined for almost every uj by l|3.22|l . 
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expressing the solution to ()3.16|) as a superposition of eigenfunctions. (There is a 
similar formula for the advanced solution with Im^ < 0, involving an upper bound 
ti on the support of the external force.) 

Now let us fix w and suppose that we have an eigenfunction of the form (|3.22p . 
Then by (|3T8)) and ((3T9l) . 

(3.25) Vc.(^,s) = '-^^ r sin{u\s'~s\)<,{x,s')ds' + ^4r-cos{uJs)^x,uj) 

(3.26) -ydlcjj^ix) = -lj'^{1 + Rexix,Lu))(j)^{x) ~ iLu'^lmxix,Lu)g^{x), 
with 

(3.27) gLo{x) = lini < <j){x , u! + iS) + (f>{x , uj — iS) > . 

6 — >0 L J 

In other words, the eigenfunctions which appear in the expansion (j3.24p are of the 
form p.9p with b{x) = 0. We refer to such eigenfunctions as spectral eigenfunctions. 
By (|3.24p , the spectral eigenfunctions form a complete set for the description of the 
dynamics of the extended system (|1.36p . Thus, the freedom to choose g^j in p.26p 
provides us with a rich enough family of solutions to describe the dynamics of the 
TDD string. 

Note that for C = w G M, a solution to (|3.16p satisfies the eigenmode equation 
away from the spatial support of the driving force. Thus, one may try to produce 
eigenfunctions via a solution to (|3.16p with a "source at infinity," that is as a limit 
of a sequence of solutions to (|3.16p with driving forces supported farther and farther 
from the origin. Depending on the susceptibility and the sequence of driving forces, 
the limit may or may exist and be a non zero. However, when this procedure leads 
to a non-trivial limit, the resulting eigenmodes are of the form (j3.26p with a special 
choice of the arbitrary function g^^. If we solve (I3.16P in the upper or lower half 
planes, we obtain in this way 

(1) The causal eigenfunctions with gu{x) = (f>i^{x), 

(3.28) ^P^ix,s) = ^^ r e'^\''~'k{s',x)ds' 

2 J-oo 

(3.29) idlM^) = -u;\l + x{x,u;))M^), 

corresponding to a driving force in the distant past and the causal boundary 
condition with the strings at rest at i = — oo. 

(2) The anti-causal eigenfunctions with guj{x) = —(j)^{x), 

(3.30) Mx,s) = ^ re--l^'-^l^(s',x)ds' 



(3.31) idlM^) - -Lu^{l + xix,co))Mx), 

corresponding to a driving force in the distant future and the anti-causal, 
or advanced, boundary condition with the strings at rest at t — 4-oo. 
Note that the causal eigenfunction equation p.29p is simply the formal eigen- 
value problem obtained from the time Fourier transform of the TDD string equation 
(|1.4p . Similarly, the anti-causal eigenfunction equation p.3ip is the formal eigen- 
value problem obtained from the time reversal of (jl.4p . However, the solutions to 
p.29p and (|3.3ip are very special in nature and need not provide a complete set for 
the expansion of solutions to (|1.4p . Indeed, it can happen that there are no physi- 
cally relevant causal or anti-causal eigenfunctions. For instance, in a homogeneous 
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string the only solutions to ()3.29p and ( 3.3ip are exponentially growing as 



X 



or X — > — oo, and are thus irrelevant^ Nonetheless, there are of course plenty of 
spectral eigenfunctions coming from the extended system. In section [S] below we 
construct a rich family of bounded plane wave solutions for a homogeneous string, 
using (|3.26p with suitable choices for the source term g^. 

This is not to say that the causal and anti-causal eigenfunctions are always ir- 
relevant. That depends on the physics, and indeed the formal eigenvalue problem 
can be useful in the right context. For instance, if dissipation is restricted to a 
proper subset of the string, or more generally if hnx^x, co) falls off sufficiently fast 
as a: — > ±cxd, then among the causal eigenfunctions are scattering solutions which 
describe the reflection and transmission of a plane wave emitted from a source at 
X — ±cxD. In the complete family of spectral eigenfunctions these solutions are very 
special, however they are precisely those solutions needed to analyze the extended 
dynamics (|1.36p for a driving force situated very far from the dissipative portion 
of the string. In section [6] we illustrate the role of causal eigenfunctions in scatter- 
ing theory by computing the scattering modes for a string in which dissipation is 
restricted to a; > 0. 

4. Energy flux in an eigenfunction 

The energy density £ in an eigenfunction, at a point x with f(x, lu) ^ 0, is typ- 
ically infinite due to the contribution from the hidden strings. This result is to 
be expected physically as the eigenfunction represents the idealized steady propa- 
gation of a monochromatic wave through an absorbing medium with infinite heat 
capacity. However, due to the decoupling between the hidden strings, energy can 
flow only through the physical string and we expect the energy flux J to be finite. 
Furthermore as we shall see d^J, which by (|2.7p is formally equal to —dt£, can be 
non-zero at a point x, in which case the eigenfunction is a steady state in which 
energy is dissipated to or absorbed from the medium at a; at a constant rate. 

By p.Sp . the energy flux for an eigenfunction is 

(4.1) J{x,t) = -jReliLU(j)uj{x)dx<puj{x)> = 'yujlm(j)^{x)dx4>ui{x). 

Thus, by ((3:261) . 

—dxJ{x) = —ijjlm(l)^{x)^dl.4>^{x) 
(4.2) 



LoHm x{x, cj)Re U)^ {x)g^ {x) | 



for a spectral eigenfunction, where guj{x) is the arbitrary function describing the 
excitation of the medium, as represented by the hidden strings. The energy density 
of the physical string 

(4.3) Eoix) - ^{u;^\Mx)\^+j\dxCJ,Ux)\'}, 

is constant in time, and there is a constant rate of dissipation at each x with 
-dxJix) ^ 0, 

(4.4) dt£{x) = dtEhsix) - w3lm£(2;,c^)Re |^4^5^(x)}. 

■'One can use the abstract theory of Gelfand triples and generalized eigenfunction expansions to 
show that only those eigenfunctions which grow at infinity slower than (1 + |a'|)^+ with arbitrary 
<5 are relevant in the eigenvalue expansion of 1 11.3611 . 
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Thus the cigcnfunction represents a steady state situation in which energy is flowing 
into or out of the hidden strings at a constant rate at each point x with Im xC^^i ^) 7^ 
0. 

For a general spectral eigenfunction (|3.26p . the dissipation dt£{x) may be positive 
or negative, however for a causal eigenfunction (j3.29p . with g^^ — 4>u}, we have 

(4.5) dt£{x) — dtEfis{x) — uj^Iinx{x,uj)\(j)uj{x)\'^ > 0, (causal eigenfunction) 

corresponding to a steady dissipation of energy from the physical string into the 
medium, represented by the hidden strings. Similarly, for an anti-causal eigenfunc- 
tion (|3.3ip there is a steady flux of energy out of the medium and into the physical 
string 

(4.6) dtEiis{x) = — a;''lmx(a;,cj)|0i^(a;)p < 0, (anti-causal eigenfunction) . 

5. Plane waves and momentum flux in an homogeneous system 

Suppose that the susceptibility x(x, t) = xC''") is constant for the whole range of 
a; e M. Then it is interesting to look for plane wave eigenfunctions e'(''^~"*'((/)o, ttq, 
ipois),Oo{s)) with ^Oj^'o constants and ■00 (s), 6*0(5) independent of x. 

A first observation is that there are no causal or anti-causal eigenfunctions of 
this form, at least at frequencies with Imx(w) > 0. Indeed a causal solution would 
satisfy 

(5.1) -7fc''/'o = -uj^{l+xH)cj)o, 

and the only non-trivial solutions to this equation are exponentially growing as 
X —> ±cx) unless x(^) is real. Such evanescent waves play a role in constructing 
scattering states for inhomogeneous systems but are physically irrelevant in the 
homogeneous system. 

Thus to find plane wave solutions it is necessary to look beyond causal eigen- 
functions. For a plane wave spectral eigenfunction, see (I3.26p . we have 

(5.2) -lk^(l)o = -iw^5oInix('^) -^^(1 + Rex(w))0o, 

(5.3) Vo(s) = -^ cos{los)'^{uj) + -^ / sin(Lj|s' - s|).?(s') ds', 

with 50 a-ii arbitrary constant describing the excitation of the hidden strings. A 
bounded solution results only for k real, so 

(5.4) 50 = ia'/'o, a G M, 

must be a pure imaginary multiple of 0o. Furthermore, we must have 

(5.5) 1 -I- Rex(t^) - almxiuj) > 0. 
Then setting 

(5.6) 7fc2 = u;2(l-|-Rex(w)-aImx(w)), 

we obtain a non-trivial plane wave solution. 

If the medium is not dissipative at frequency uj, so Imx(w) = 0, there is a 
dispersion relation between k and w 

(5.7) 7fc2 = a;2(l + Rex(c^)), (Inix(a>) = 0). 
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At frequencies to with non-trivial dissipation, that is Ivaxi'^) 7^ 0, there is no relation 
between k and lo. Indeed k may be chosen arbitrarily provided we take 

(5 8) a = 1 + RcxM _ 7fc' 

Imx(cj) uP-\va.x(y>^ 

Observe that (|5.7p is essentially the Fourier transform in time and space of 
the formal eigenvalue problem (|5.ip in the case that Imx(a;) = 0. The lack of a 
dispersion relation between k and a; in a dissipative medium highlights the failure of 
the formal eigenvalue problem in this context. Indeed, we have seen that the formal 
eigenvalue problem (|5.ip has no physically relevant solutions in a homogeneous 
medium. Nonetheless, there is a complete eigenvalue basis for the extended system 
(|1.36p consisting of the plane wave solutions derived above. 

By (|4.ip . the energy flux for a plane wave 

(5.9) J = 7^fc|'/'o|' 

is constant. Thus dxJ — and plane wave solutions represent a steady state 
flow without dissipation of energy. More precisely, the dissipation of energy to 
the medium, as described by the hidden strings, is exactly balanced by the energy 
emitted from the medium. The energy density of the physical string is 

(5.10) Ea = {lu^ + -/k^) l^ol'' 

and the energy density of the medium, described by the hidden strings, is of course 
infinite. 

As the system is homogeneous, we can also consider the wave momentum density 
and stress in a plane wave eigenfunction. By (|2.17p the wave momentum density 
of the physical string is 



Rexito)). 

UJ 

The wave momentum density of the hidden strings is infinite, since 
(5.14) phs = iok j \Ms)\''As 



(5.11) 




Po = wfc ( 


■poY -fcImAo(/)o, 


(5.12) 


Ao = 


= / <r(s)-0o(s)ds : 


= iw(/)o (alnix(w) 


Thus 








(5.13) 


Pq = 


= cjfc{l + Rex(w)- 


-almx(w)}|0|2 : 



by dUHD. 

Of more interest is the stress, which by (|2.191 I2.20p is 

(5.15) T = Ea+Lolu{Ka(j}a+Ths. 
where Eq given by (|5.10p and Ths is formally 

1 f°° 

(5.16) Ths ^ T^ {c^Vo(5)|'-|9.Vo(s)|2}ds. 

The integral on the r.h.s. of (|5.16p does not converge absolutely, but we can regu- 
larize it by defining Ths — linr^^o Ths{5) with 

1 1'°° 

(5.17) Ths{S) = - e~'\^^ {l,'\Ms)\^ - \dsMs)\''} ds- 
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To evaluate the integral in (|5.17|) it is useful to write 

(5.18) 2 

where we have used p. 31) . Thus 
(5.19) 

= l\MO)\'-^j e-^l^i|Vo(s)pds--c.y e-'\^k{s)ImMs)cj,ods. 

It follows that 

1 f°° 1 

(5.20) Ths = limr,is((5) = --tj / <r(s)ImV'o(s)0ods = --tjlm Ao(^o, 

5^0 A J_oo A 

if, for instance, \'^{s)\ is integrable. Thus 

T == £'o + -cjImAo^o 

(5.21) 2 

which is identical to the result for the undamped string! 

6. Scattering eigenfunctions for a semi-infinite medium 

To close, we would like to illustrate the power of the above method by unam- 
biguously constructing the scattering states representing an idealized description of 
the following experiment. Imagine that we have an long, i.e., infinite, string whose 
right half (x > 0) is subject to dispersion and dissipation with susceptibility %('''). 
That is, the the susceptibility for the whole string is 

(6.1) ^(-'-)-|"r^ '^o- 

Suppose we drive the left end of the string, at x = — oo, periodically with frequency 
w, sending an incoming wave toward the TDD right half of the string. After some 
time a steady state is reached in which a certain fraction of the incoming wave is 
absorbed by the dissipative half of the string and a certain fraction is reflected. 

The steady state eigenfunction describing the above experiment is a causal eigen- 
function since the source is at a; = -co and in the distant past. Thus by (j3.29p . 
the string displacement satisfies 

(6.2) ldlK{^) = -c^2(l + /[x>O]x(c^))0^(a;), 

where I\x > 0] = 1 for x > and otherwise. Indeed, this is the naive equation 
that one might right down for the scattering states of a TDD string. However, it 
is only within the context of the extended Hamiltonian system that we understand 
this to be only one of many possible eigenfunction equations, the choice of which is 
dictated by the physics under consideration, namely a source at spatial infinity in 
the distant past. 
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By (|6.2|) . we have 

(6.3) 0^(x) = e"=<" + re-"=<% fc< = ^, {x < 0) 

on the left half of the string, up to an over all multiple which we fix to be 1 without 
loss. In (|6.3p the term e''^<^ is the incoming wave and re~''^>^ is the reflected wave. 
To find the reflection coefficient r we need to solve for the form of the eigenfunc- 
tion on the right half of the string. Again by (|6.2I) . we have 

(6.4) ldl.pUx) = -uj\l + x{to))Mx), (a:>0) 

for a; > 0. Furthermore, since we require the solution to be bounded this determines 
(p uniquely up to an over all multiple 

(6.5) cj,^^) = ^;e"=>^ 

where v is an, as yet, undetermined transmission coefficient and 



(6.6) fc> = +J^(l + x{u;)). 

Here ^{/z denotes the unique square root of z in the upper half planey Since 
Imfc> > 0, at least if Imx(a;) ^ 0, we see that 4>u}{x) decays exponentially in the 
dissipative half of the physical string (x > 0) . The excitation of the hidden strings, 
which are restricted to a; > 0, is given by p.28p : 

e-k-«\(s')ds', (x>0). 

-oo 

To determine v and r we note that the eigenfunction equation (|6.2p forces (f)^ 
and dx4>uj to be continuous functions of x, in particular at a; = 0. Thus, 

(6.8) lim(/)„(a;) — 1 ^ r = v = \vca(l)^{x) 

(6.9) \\mdx(t>i^{x) = fc<(l-r) = kyv = \\mdx4>i^{x). 

x]0 xiO 

We conclude that 

(6.10) V 
where 



2 _ l-/3(w) 



(6.11) 




yi + xH w>o, 

- v^TxH t^ < 0. 



It is useful to consider some general properties of p{lo) and the reflection and 
transmission coefficients. To begin note that since toluYxi'^) > we have 



(6.12) ujRe yi + x(w) > 0, 



If x{^) is I'^al and (1 + xC'^)) > then the medium transmits at frequency uj and A;> should 
be determined as the Umit 



fc> = limt/^^^^i^(l + x(^ + i5)). 
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and thus 

(6.13) Rep{uj) > 0. 
This iniphes that 

4 Rpfr^ 

(6.14) l-|r|2 = ,., Rep{cj) = \v\'Rep{u;) = \vf^ > 0, 

so, in particular, 

(6.15) |r| < 1. 

Eq. (|6.14p expresses the continuity of the energy flux J(x) at a; = 0, since by 
(SU) 



J{x) = ^Lulm(J3i^{x)dx4>uj{x) 

ffc<(l-|r|2) x<0 

(6.16) ~ ^'^ JRe/c>|v|2e-2i'"'=>^ x>0 

Furthermore, we see that the energy flux is non-negative, representing a flow of 
energy from the source at a; = — oo, and the rate of dissipation, —dxJ{x) ^ 0, is 
non-zero for every x > 0, 

(6.17) -dxJ{x) = 2^7^2(1 -|rnim/c>e-2inifc>^, (x > 0). 
Finally, we note that the energy density of the physical string is 

Eo{x) = ]^{uj^\cj>^{x)\'' +-i\dx<i>^{x)\^) 

= u\l + \r\^)[^ ,^ , '^<°' 

where we have noted that by (|6.10p 

(6.19) 1 + l.p = ^ + '^^("^'^.1^ = i±E±lMH^. 



In summary, the scattering eigenmode 

e'"(;^-*) + re""^V^-V x <Q 



-t] , _-i^(-^+t) 



(6.20) c^{x,t) = , . 

ue V v^ ; a; > 

(6.21) ^{x,s,t) = ^e'"('^^*) /e'"l^-^\(s')ds', (a; > 0), 

describes a steady state in which the incoming wave is partially reflected, with the 
remainder an evanescent transmitted wave that penetrates the dissipative part of 
the string with an exponential profile resulting in an excitation of the hidden strings 
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accounting for dispersion and dissipation. The total rate of dissipation in the TDD 
portion of the string is 

/•CO 

(6.22) - / d^J{x)dx = J(0) - V7^'(l - |r|2). 

Jo 
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